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We study analytically the corrections to the leading terms in the Renyi entropy of a massive 
lattice theory, showing significant deviations from naive expectations. In particular, we show that 
finite size and finite mass effects give rise to different contributions (with different exponents) and 
thus violate a simple scaling argument. 

In the specific, we look at the entanglement entropy of a bipartite XYZ spin-1/2 chain in its 
ground state. When the system is divided into two semi-infinite half-chains, we have an analytical 
expression of the Renyi entropy as a function of a single mass parameter. In the scaling limit, we 
show that the entropy as a function of the correlation length formally coincides with that of a bulk 
Ising model. This should be compared with the fact that, at criticality, the model is described by 
a c = 1 conformal field theory and the corrections to the entropy due to finite size effects show 
exponents depending on the compactification radius of the theory. We will argue that there is no 
contradiction between these statements. 

If the lattice spacing is retained finite, the relation between the mass parameter and the correlation 
length generates new subleading terms in the entropy, whose form is path-dependent in phase-space 
and whose interpretation within a field theory is not available yet. These contributions arise as 
a consequence of the existence of stable bound states and are thus a distinctive feature of truly 
interacting theories, such as the XYZ chain. 
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I. INTRODUCTION 

In the past few years, there has been a growing in- 
terest in quantifying the degree of "quantum-ness" of a 
many-body state, which is usually taken as the ground 
state |0) of a given HamiltonianirJ,. As the entangle- 
ment constitutes an intrinsically quantum property, one 
popular way to measure this is through bipartite entan- 
glement entropy^*. To this end, the system is divided 
into two subsystems (A and B) and one looks at the re- 
duced density matrix, obtained by tracing out one of the 
two subsystems 

p A = Tr B \Q)(0\ . (1) 

We will be interested in the Renyi entropies^ 
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lnTr,32 



(2) 



which, in the a — > 1 limit, reduces to the Von Neumann 
entropy 



-Trp^log p A 



(3) 



Varying the parameter a in @ gives us access to a lot of 
information on pa, including its full spectrum^. 

For gapped systems, the entanglement entropy satisfies 
the so-called area law, which means that its leading con- 
tribution for sufficiently large subsystems is proportional 
to the area of the boundary separating system A from B. 



In 1 + 1 dimensional systems, the area law implies that 
the entropy asymptotically saturates to a constant (the 
boundary between regions being made just by isolated 
points). 

Critical systems can present deviations from the simple 
area law. In one dimension, in particular, the entangle- 
ment entropy of systems in the universality class of a 
conformal field theory (CFT) is known to diverge loga- 
rithmically with the subsystem size^ilii. From CFT, a lot 
is known also about the subleading corrections, which, in 
general, take the unusual formii~— 
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1 + - )\n — + c' a + b a (£) £~ 2h/a + 
a J a 



(4) 

where c is the central charge of the CFT, £ is the length 
of subsystem A, ao is a short distance cutoff, c' a and 
b a (£) are non- universal and the latter includes a peri- 
odic function of £ (with the period given by the Fermi 
momentu m 1 2 ' 1 4 ) , and h is the scaling dimension of the 
operator responsible for the correction (relevant or irrel- 
evant, but not marginal, since these operators generate a 
different kind of correction, which will be discussed later) . 
This result is achieved using replicas, and thus, strictly 
speaking, requires a to be an integer. Moreover, it should 
be noted that in Ref.— the corrections are obtained from 
dimensionality arguments, by regularizing divergent cor- 
relations by an ultra-violet cut-off ao- Thus, technically, 
the subleading contributions in Q are extracted from 
scaling properties and are all of the form £ja§. 
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Determining the exponents of the corrections is im- 
portant both in fitting numerics (where often really large 
£ are unobtainable) and also for a better understanding 
of the model. For instance, the scaling exponent h also 
determines the large n limit of the entropy (single copy 
entanglement)^. Moreover, especially for c = 1 theories, 
h provides a measure of the compactification radius of 
the theory^ and thus of the decaying of the correlation 
functions. Up to now, this conjecture has been checked 
in a variety of critical quantum spin chains models^—. 

Moving away from a conformal point, in the gapped 
phase universality still holds for sufficiently small rele- 
vant perturbations. Simple scaling arguments guarantee 
that the leading terms survive, but with the correlation 
length £ replacing the infra-red length-scale I. Recent 
results, based on exactly solvable models, indicate the 
appearance of the same kind of unusual corrections to 
the Renyi entropies, which are now functions of the cor- 
relation length £ with the same exponent h^: 



S a =^(^)ln^+A a + B a r h / a 
12 V a J a 



(5) 



There is a factor of 2 difference in each term between (|4]) 
and jl]), due to the fact that the first is a bulk theory 
(with both chiralities in the CFT), while the latter is 
expected to be akin to a boundary theory, were only one 
chirality in the light-cone modes effectively survives. 

In this paper, we are going to investigate the sublead- 
ing terms in the Renyi entropy of the one-dimensional 
XYZ model. In the scaling limit we find that the entropy 
is (modulo a multiplicative redefinition of the correlation 
length) 
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where (7_i(n) is a divisor function, defined by (|41[) . 

While the leading term correctly reproduces a c = 1 
central charge, the interpretation of the scaling exponents 
in the subleading addenda is less straightforward. Com- 
paring the exponent of the first correction with ([5]) would 
indicate h = 2, i.e. a marginal operator. This term 
cannot be due to the same operator acting at the crit- 
ical point since, as shown in RefJ^, marginal operators 
give rise to logarithmic corrections. Moreover, the critical 
XXZ chain is known to have relevant fields with h = K 
(the Luttinger parameter) and the opening of a gap in 
the XYZ model implies the presence of a h < 2 operator. 
We will show that the operator content that can be ex- 
tracted from ([6]) matches that of a bulk Ising model and 
the first correction can be interpreted as arising from the 
energy field. Notice that the leading, logarithmic term 
can thus be equally interpreted as with c = c = i . 



Furthermore, if we include also lattice effects, which 
vanish in the strict scaling limit, additional corrections 
appear in ([6]) and, while they are less important than 
the dominant one for sufficiently large a, they can be 
relevant for numerical simulations in certain ranges of a. 
These corrections turn out to be path dependent (prob- 
ably due to the action of different operators) and many 
kind of terms can arise, such as £ -2 ' 1 , ^~ 2h / a , £-( 2 -' l ) 5 
£-2(1+1/0!) or cven l/ln(£). In light of RefJ^, some of 
these terms were to be expected, but the others still lack 
a field theoretical interpretation, which might be possi- 
ble by applying a reasoning similar to that of RefJ^ for 
a sine-Gordon model. 

The paper is organized as follows: in Sec. [II] we intro- 
duce the XYZ model, its phase-diagram and elementary 
excitations. In Sec. Mil we discuss the structure of the 
reduced density matrix obtained by tracing out of the 
ground state half of the system and its formal equiva- 
lence with the characters of an Ising model. In Sec. IIVI 
we present the full expansion of the Renyi entropy and 
discuss how in the scaling limit this expansion coincides 
with that of a bulk Ising model. Then we discuss the lat- 
tice corrections both in the regime where the low-energy 
excitations are free particles and where they are bound 
states. In the latter case, we show that different paths of 
approach to the critical point give rise to different sub- 
subleading corrections. Finally, in Sec. |V]we discuss our 
results and their meaning. 



II. THE XYZ SPIN CHAIN 

We will consider the quantum spin-i ferromagnetic 
XYZ chain, which is described by the following Hamilto- 
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(7) 

where the a" (a = x, y, z) are the Pauli matrices acting 
on the site n and the sum ranges over all sites n of the 
chain; the constant J, which we take to be positive, is an 
overall energy scale while (J y ,J z ) take into account the 
degree of anisotropy of the model. 

In Fig. [1] we draw a cartoon of the phase diagram of 
the XYZ chain. The model is symmetric under reflections 
along the diagonals in the (J Zl J y ) plane. The system is 
gapped in the whole plane, except for six critical half- 
lines/segments: J z = ±1, \J y \ < 1; J y = ±1, \J Z \ < 1 
and J z = ±J y , \J Z \ > 1. All of these lines correspond 
to the paramagnetic phase of an XXZ chain, but with 
the anisotropy along different directions. Thus, in the 
scaling limit they are described by a c = 1 CFT, with 
compactification radius varying along the line. We will 
use < /3 < V87r (the sine- Gordon parameter of the cor- 
responding massive theory) to parametrize the radius. 
The critical segments meet three by three at four "tri- 
critical" points. At these endpoints, f3 is the same along 
each line, but the different phases have a rotated order 
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FIG. 1: (Colored online) Phase Diagram of the XYZ model 
in the (J z , J y ) plane. The blue solid lines -J x = ±1, \ J y \ < 1; 
J y = ±1, \J,\ < 1 and J z — ±J y , \J X \ > 1- correspond to 
the critical phase of a rotated XXZ chain. Out of the four 
"tricritical" points, Ci,2 are conformal and Ei,2 are not. The 
area within the red rectangle is the portion of the phase dia- 
gram we study in this article, and it is completely equivalent 
to the principal regime of the eight-vertex model. The yellow 
and green lines are the curves of constant I and n respectively 
in the (J z , J v ) plane, according to the parametrization {TTJ. 



parameter. Two of these points -C\ % % = (1, —1), (—1, 1)- 
are conformal points with /3 2 = 87r; while the other two 
-Ei t 2 = (1, 1), (— 1, — 1)- correspond to /3 2 = and are 
nonconformal. The former points correspond to an anti- 
ferromagnetic (AFM) Heisenberg chain at the BKT tran- 
sition, while the latter describe the Heisenberg ferromag- 
net. Thus, at E\2 the system undergoes a transition in 
which the ground state passes from a disordered state to 
a fully aligned one. Exactly at the transition, the ground 
state is highly degenerate while the low-energy excita- 
tions arc gaplcss magnons with a quadratic dispersion 
relation. 

In studying the XYZ chain, one takes advantage of 
the fact that (J7J commutes with the transfer matrices 
of the the zero-field eight- vertex model (see for example 
Refs i 18 i 19 ) and thus the two systems can be solved simul- 
taneously. The solution of the latter is achieved through 
the parametrization of J y , J z in terms of elliptic functions 
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cn(iA) dn(iA) 
l-k sn 2 (iA) 
1 + fc sn 2 (iA) 
1-fc sn 2 (iA) 



(8) 



where (r, A) are the well-known Baxter parameters^, 
and sn(x), cn(x) and dn(x) are Jacobian elliptic functions 
of parameter k. A and k are parameters, whose natural 



domains arc 



< k < 1 



< A < I(k') 



(9) 



J(fc') being the complete elliptic integral of the first kind 
of argument k' = \/l — k 2 . 

The definition of (A, T) itself is particularly suitable to 
describe the anti-ferroelectric phase of the eight-vertex 
model (also referred to as the principal regime), corre- 
sponding to A < —1 and |T| < 1. However, using the 
symmetries of the model and the freedom under the re- 
arrangement of parameters, it can be applied to the whole 
of the phase diagram of the spin Hamiltonian (for more 
details see RefJ^). For the sake of simplicity, in this pa- 
per we will focus only on the rotated principal regime: 
Jy > 1, \Jz\ < 1, sec Fig. [TJ and we defer to a different 
publication some interesting properties of the generaliza- 
tion to the whole phase diagram. 

Before we proceed, it is more convenient to switch to 
an elliptic parametrization equivalent to ([8]) : 
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The elliptic parameter I corresponds to a gnome r 
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which is half of the original. The rela- 
tion between k and I is known as Landen transformation. 
Note that < \i < ir. 

In terms of these new parameters, we have 
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dn[2iI(I')A0r; Z] ' 



cn[2H(l')fi/Tv;l] 
~ dn[2il(l')n/ir;l] 



(11) 

Curves of constant / always run from the AFM Heisen- 
berg point at fx = to the isotropic ferromagnetic point 
at fi = ir. For I = 1 the curve coincides with one of the 
critical lines discussed above, while for / = the curve 
run away from the critical one to infinity and then back. 
In Fig. [T] we draw these curves for some values of the 
parameters. 

For later convenience, we also introduce 



x = exp 
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27(A) 



= e' 1 



(12) 



Using Jacobi's theta functions, we introduce the elliptic 
parameter k± connected to x, i.e. 



h = 



j(0, *) x?(-l;x 2 ) 
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k(x) , (13) 



or, equivalently, ^ ™ — i^ 7 " (i- c -i ^ is to r what k\ is 
to ht/tt). In (fTB"]) we also used the (j^Pochhammer symbol 



n-l 



(a;q)„ = ]J(l-aq k ) 



(14) 



The correlation length and the low-energy excitations 
of the XYZ chain were calculated in Refill. There are 
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two types of excitations. The first can be characterized 
as free quasi-particles (spinons). The lowest band is a 
2-paramctcr continuum with 



AE llcc (q 1 ,q 2 ) 



1 — k( cos 2 qi + d 1 — k 2 cos 2 q 2 



The energy minimum of these state is achieved for (71,2 
0, ±7r and gives a mass gap 
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For /j, > 7r/2, in addition to the free states just dis- 
cussed, some bound states become progressively stable. 
They are characterized by the following dispersion rela- 
tion 



AE s (q) = -2J- 



[21(1')^; I'} 1(h) 
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sn(sy;fci) 



1 — dn 2 (sy; fc^) cos 2 — 



1 — cn 2 (sy; fe^) cos 2 — , 



(17) 



where y = il(ki)r (^ — l) and s counts the number of 
quasi-momenta in the string state. In the scaling limit, 
these bound states become breathers. The mass-gap for 
the bound states is (setting q = above) 



AE S = AE b00 sn(sy;k[) , 



(18) 



from which one sees that for (j, > tt/2 the s = 1 bound 
state becomes the lightest excitation. 

The correlation length for the XYZ chain (Fig. ^ was 
also calculated in Reh2& and it is given by: 



■\ mfc 2 

2 dn 2 [i2/(fe 2 )^( M -f );fe^j 



< (J, < f . 

f < LI < 7T 



(19) 

where ao is a short distance cut-off, such as the lattice 
spacing, that sets the length unit and the new parameter 
k 2 is the Landen transformed of k\ : 



k 2 = k(x 2 ) = 



1 - k{ 
1 + k'i 



(20) 




FIG. 2: Curves of constant correlation length of the XYZ 
model in the (J z , J y ) plane. Regions of similar colors corre- 
spond to the same correlation length values and red colors are 
associated to higher values of £. It is worth stressing that the 
correlation length does not show any essential critical behav- 
ior, whilst the block entropy shows it in proximity of points 
Ei and E2. 



limit where the infinite chain is partitioned in two (semi- 
infinite) half lines. For this configuration, the reduced 
density matrix can be computed as the product of the 
four corner transfer matrices (CTM) of the correspond- 
ing eight-vertex model 2 ^— . In Ref.— it was shown that 
it can be written as 



1 

x 2j 



(21) 



where Z = (—x 2 1 x 2 ) 00 , that is, the partition function of 
the eight-vertex model, is the normalization factor that 
ensures that Trp — 1. Thus we have 



Trp a = 



(~x 2 ]x 2 ) a x 



(22) 



The first behavior in flT9l is due to the free particles 
states, while the s — 1 bound state is responsible for the 
second 2 ^. 



III. THE REDUCED DENSITY MATRIX 



The bipartite entanglement entropy for the ground 
state of the XYZ chain was calculated i n 22 i 23 , in the 



and for the Renyi entropy 



S a — 
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00 
3=1 



ln(l 



The structure (Eqs. [2U 
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EM 1 



Ma \ 



(23) 

for the reduced density 



matrix of the half-line is common to all intcgrable, local 
spin- 1/2 chains^ and thus the entanglement spectrum of 
these models is the same and only depends on x, which 
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FIG. 3: Curves of constant entropy of the XYZ model in 
the (J z ,Jy) plane. Regions of similar colors have similar en- 
tropy values and the line where colors change are the lines 
of constant entropy. The brighter is the color, the bigger the 
entropy. 



in this context is usually parametrized as x = e~ e . For 
the XYZ chain, e = — i/ir. In Rcf.— we showed that 
e (and thus the entropy) has an essential singularity at 
non-conformal points, and thus its behavior differs dra- 
matically from the conformal one. 

From (|2"3")h one can see that the Renyi entropy is a 
monotonically decreasing function of e: 



Using ||T2j} and iJ 



lim S a = oo , 
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lim S a = . 
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(24) 



(25) 



we can plot the entanglement entropy in the phase di- 
agram of the XYZ model. In Fig [3] we show a contour 
plot of the Von Neumann entropy in the (J z , J y ) plane, 
from which one can clearly see the different behavior of 
the conformal and non-conformal points. 

g-products of the form (f22j) give easy access to the 
spectral distribution of the reduced density matrix, since 



fc=l 



q 



= 1 



2^p 

n=l 



(1) 



(26) 



where p^(n) is the number of partitions of n in distinct 
positive integers. Also note that, since 



k=l 



IK 1 

k=l 



„2fc-l\-l 



(27) 



p^(n) — po(n), that is, the number of partitions of n 
into positive odd integers. 

Moreover, one can recognize Z to be formally equal to 
the character of the spin field of the Ising CFT. To show 
this, we write 

oo oo - 2j 

2(^x 2 )=n( i +^)=riT^r ( 2s ) 

3=1 i=i 1 

and we use the Euler's formula for pentagonal num- 
bers (which is a consequence of the Jacobi triple-product 
Identity)^ 



to recognize that 
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Z = X-^^ (k/TT) 



(30) 



where Xi S 2 nS ( T ) i s the character of the spin = 1/16 
operator of a c = 1/2 CFT: 



X^f'V) = 
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model. We have 



with q = c ImT and (p,p') = (4,3) for the Ising minimal 



Trp Q 



Ising/- / \ 



(32) 



As we discussed above, the critical line of the XXZ 
chain is approached for I — > 1, that is, for r — > and 
;r — > 1. On this line, excitations are gapless and in the 
scaling limit the theory can be described by a conformal 
field theory with central charge c = 1 . To each < fj, < ir 
it correspond a different point on the critical line, with 
sine-Gordon parameter 1 = 8n (l — H\28, The two end- 
points are exceptions, since fj, = 0, tt identify the same 
points for every /. However, while in the conformal one 
we still have x — > 1, close to the ferromagnetic point, 
around /.t = 7r, both x and q can take any value between 
and 1. Hence a very different behavior of the entangle- 
ment entropy follows ! 23 ' 29 

In order to study the asymptotic behavior of the en- 
tanglement entropy close to the conformal points, it is 
convenient to use the dual variable 



(33) 



which is such that x — >• as I — >• 1. 
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Expressions like (|22j) involving q-products can be writ- 
ten in terms of elliptic theta functions. For instance, 
this was done for the entanglement entropy of the XYZ 
To study the conformal limit, one performs a 



,23,30 



modular transformation that switches 



x. Since 



k(x) = fc'(x) = k[ = 
and using (|13p . we have 

(-* 2a ;* 2a L = 

1/12 



(x,x ) 00 



k 2 {x a ) 



(34) 
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(36) 



The modular transformation that allowed us to switch 
from x to x is the same one that connects characters 
in minimal model of inverse temperature. For the spin 
operator of the Ising model we have 



Ising / \ 
Xl,2 ( T ) 
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Ising/ 1 / \ Ising/ n / \ 

Xi,i (-V r )-X2,i (- 1 /' r ) 



(37) 



Using (|3Tj) and the identities (|27l and [29]) one can prove 
that 



IV. EXPANSIONS OF THE ENTANGLEMENT 
ENTROPY 

Close to the conformal points, Eq. (f2"3")l is just a formal 
series, since x ~ 1. However, using (p?6"|) . it is straight- 
forward to write a series expansion for the Renyi entropy 
@ in powers of x <C 1: 
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where the coefficients 
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= - E (i+ fc )+ E - 



cri(n) 



(41) 
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is a divisor function^ 3 , and takes into account the expan- 
sion of the logarithm over a q-product and play a role 
similar to the partitions of integers in (121)1) . It is worth 



noticing that the constant term ln(2 -1 / 2 ) = \n(S^ 16 ) - 
where S^ le is an element of the modular S-matrix of the 
Ising model - is the contribution to the entropy due to 
the boundary 3 ^. 

The a — > 1 yields the Von Neumann entropy: 



S a = - — In x - - In 2 



(42) 
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which agrees with ([36ft and implies 



Ising /• 

=1 Xi,i s (i7r/ae) ~ X 2 ,i {w/ae) 



Trp a = 2— 



Ising/ 



Ising/- 

Xi.i (i7r/ e )-X 2 .i (w/e) 



Ising 



(38) We see that, contrary to what happens for a > 1, all sub- 
leading terms -which are powers of x- acquire a logarith- 
mic correction, which strictly vanishes only at the criti- 
cal points. A more specific analysis of this phenomenon, 
together with a quantitative study of what introduced 

^3g^ in the next subsections, will be the subject of a future, 
partly numerical, publication. 



This agrees with what conjectured in RefJ^, but with 
the important difference that the characters in (|39|) arc 
c = 1/2 and do not belong to the infrared c = 1 bulk de- 
scription of the XYZ chain^i. This Ising character struc- 
ture for the CTM of the eight-vertex model was already 
noticed, see Refi^. We acknowledge that, once the for- 
mal equivalence (fBT))) is understood, the content of the 
last page follows almost trivially, but we decided to pro- 
vide a brief derivation here for the sake of completeness 
and because this result has consequences on the struc- 
ture of the entanglement entropy, a fact which is not well 
known and will be discussed in the next paragraph. Fi- 
nally, we notice that, being only a formal equivalence, 
cq. (f30|) docs not imply any underlying Virasoro algebra 
at work for CTM (as far as we know) and it is thus im- 
portant to recognize that these manipulations stand on 
more general mathematical concepts. 



A. Scaling limit 

Comparing with ([5|) and coherently with ([32]) , (|40| and 
(|4"2")l can be identified with the expansion of a c = 1/2 
theory. However, the parameter of this expansion, x, 
has meaning only within Baxter's paramctrization of the 
model ([5]). To gain generality, the entropy is normally 
measured as a function of a universal parameter, such as 
the correlation length or the mass gap. 

In the scaling limit, up to a multiplicative constant, 
one has: 



(43) 



so that, substituting this into (|4T)|). we get ([5]). Relation 
(|43l) is crucial in turning the leading coefficient in the 
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entropy of a c = 1/2 entropy such as (|40|) into that of 
a c = 1 theory, but it also doubles all the exponents of 
the subdominant corrections. It is reasonable to assume 
that this c = 1 model is some sort of double Ising, but its 
operator content does not seem to match any reasonable 
c = 1 model, since only even exponent states exist. More- 
over, comparing ([5]) with ([5]), one would conclude that a 
h = 2 operator is responsible for the first correction. It 
was argued in Ref that a marginal field gives rise to 
logarithmic corrections in the entropy, thus, either this 
corrections is due to descendant of the identity (namely, 
the stress-energy tensor), or we should think of it as a 
2h/a, with h = 1. 

In fact, we can write the partition function of the eight- 
vertex model as a bulk Ising model (i.e., quadratic in 
characters). Starting from (|38p . we have 



1 

72 



2k-\ 



k=l 



X 12 

~7T 



X 12 

~7T 



x\7 g [^)+xlT{^ 



-Isine 

Xi,i 



\Xo\ 



-ln£ 

7T 



(44) 



[Xl/2 | ~ XOXl/2 + X1/2X0 



This formulation provides a simple explanation of the 
Renyi entropy expansion ([5]) and its operator content. In 
fact, it interprets the first correction as the Ising energy 
operator, and not as a descendant of the identity. 
It also means that the prcfactor in front of the logarithm 
in the entropy can be interpreted as with c — c = i . 

No fundamental reason is known for which CTM spec- 
tra (and partition functions) of integrablc models can be 
written as characters in terms of the mass parameter x 
or x. This is the case also for Eq. (|4"4"|) . Thus, so far, we 
can only bring forth this observation while any connec- 
tion with some underlying Virasoro algebra remains to be 
discovered. To the contrary, sufficiently close to a critical 
point, the CTM construction can be seen as a bound- 
ary CFT and thus its character structure as function of 
the size of the system is dictated by the neighboring fix 
points 

We are led to conclude that, in the scaling limit, the 
entropy can be written as a function of two variables: 

S a (—, — )• When the inverse mass is larger that the 
subsystem size, we have the usual expansion of the form 
But when the correlation length becomes the in- 
frared cut-off scale, apparently a different expansion is 
possible, which, unlike (0, can contain terms with differ- 
ent exponents, like in ((6]). Hence, while the leading uni- 
versal behavior has always the same numerical value and 
scales like the logarithm of the relevant infra-red scale, 
the exponents of the corrections might be in principle 
different for terms in £ and in £. 

In the scaling limit, ao — > 0, J — > 00, and x — > in 
such a way to keep physical quantities finite. In this limit, 



only the scaling relation (|43|) survives. However, at any 
finite lattice spacing ao, there will be corrections which 
feed back into the entropy and that can be relevant for 
numerical simulations. To discuss these subleading terms 
we have to consider two regimes separately. 



B. Free Excitations: < fi < § 

For < n < J the lowest energy states are free, with 
dispersion relation (|16[) . To express the entropy as a 
function of the correlation length we need to invert (|19p . 
This cannot be done in closed form. So we have to first 
expand (|19[) . finding 



^ = 4£ 1 / 2 ^a_ 1 (2n + l)i» 

= 4^/2 + 16. 3/2 + 24 . 5/2 
3 5 



(45) 



and then to invert this (by hand) to the desired order: 

(46) 



J_a§ 

16 e 



1 _ I £0 + J_ 4 + U-6\ 

6 e 144 £ 4 ^ ' 



We get: 



1 



a inl 
ao 



1 - 2a 



On 



12a 

+B a r- +c a r 2 

-aB a C 2 ~ aB' a i 



In 2 



(47) 



where the coefficients only depend on a and contain the 
proper power of ao to keep each term dimcnsionlcss [for 

instance B a = (^f-) 2 ^]- We note that a new term 
has appeared (and more will be seen at higher orders) 
and that it is not of the forms discussed in RefJ^. 

It is worth noticing that if we express the mass-gap ([To]) 
as function of x, we would get a different series expan- 
sion and thus different corrections to the entropy These 
subleading terms would have a different form, compared 
to (|47p. and even be path dependent on how one ap- 
proaches the critical point. We prefer not to dwell into 
these details now, postponing the description of this kind 
of path-depended behavior with the bound state's corre- 
lation length to the next section. 



C. Bound states: ^ < /j, < n 

For fi > 7r/2, bound states become stable, and the 
lightest excitation becomes the s = 1 state with disper- 
sion relation ([T7|) . Accordingly, the expression for the 
correlation length is different in this region from before. 
Using the dual variable x and the formulation of elliptic 
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FIG. 4: Three different ways to approach the critical line: 
along the ^-constant lines (blue line), which is referred to as 
renormalization group flow in the text; along straight lines 
(red line); along lines that approach the criticality with zero 
derivative (green line). 



functions as infinite products, wc can write (|19l) as 

V / oo V / oo 



a 



In 



( x 2 ; x ) (x i ;x) 



(48) 



oo oo 1 

4VV — 

^ ^ 2k - 1 



cos 



71=1 k=l 



. (2r»-l)(2fc-l) 



The major difference in (|48| compared to (j45|) is that 
the bound state correlation length does not depend on 
x alone, but separately on /i and r. This means that in 
inverting ([48]) to find x as a function of £, we have to 
first specify a relation between /x and r, i.e. to choose a 
path of approach to the critical line. We will follow three 
different paths, that are shown in Fig. [4] 

(a) Renormalization group flow: The first natural path 
is (represented in the blue line in Fig.|4j: 



H = Mo 



(49) 



where s — > guides our approach to the gaplcss point. 
This path, keeping /i fixed, corresponds to the RG 
flow. In the scaling limit, the XYZ chain is described 
by a sine-Gordon model, where h is proportional to 
the compactification radius^. Thus, assuming (|4"9")) 
means changing the bare mass scale, without touch- 
ing (3. In the (J 2 , J y ) plane, this path asymptotically 
crosses the critical line with slope m = — 2/ cos/io- 
Substituting this in (|4"8")l we get: 

f = 4 5 (Mo) i 1/2 + f 9 3 (Ho) * 3/2 + O (i 5 / 2 ) , (50) 

where g{n) = cos 5-. Comparing ([50)) with the free 
case (|45p we immediately conclude that the entropy 



retains an expansion similar to (|47[) . with the differ- 
ence that all the coefficients now depend on ho and 
thus change along the critical line: 



1 



Q 



12a 



ln£ + 4 Q (//(,) + £a(/*o)C 



(51) 



(b) Straight lines in {J Z ,J V ) space: Let us now ap- 
proach a conformal critical point exactly linearly in 
the ( J z , Jy) plane: 



J v = 1 + m ■ s 



Jz 



COS flo 



(52) 



This path corresponds to the following parametriza- 
tion of t and fj, (an example of which is the red line 
in Fig. HJ): 



Hi) 



o< 



H = Ho + r(m, Ho) ■ s + 0{s 2 ) 



(53) 



where r(m, /i) = ^ ■ Thus, in the limit s — > 

0, the entropy parameter x vanishes like x cx s 71 ^ . 
Using O in (gHD 



— ~4<?(^ )i 2 +4r(m, A10 ).9'(Mo)5 2+ - +••■ (54) 

Inverting this relation, we arrive at the following ex- 
pansion of the Renyi entropy along (f5"2"|) 



1 + a 
12a 



\n£ + A a (no) 



(55) 



+B a (po)C 2/a +D a (m, f io)C 2tM ' / * 



We notice that the last term yields a new type of cor- 
rection, with a non-constant exponent, which varies 
with /1. This term is of the form £~( 2 ~' l ) ; where h 
here is the scaling dimension of the vertex operator 
c 1 ^ of the underlying sine-Gordon theory. To the 
best of our knowledge, this is the first time that such 
a correction in the Renyi entropy of gapped systems 
is discussed and it also differs from those discussed 
in RefJ^. It is surprising to see the appearance of 
the operator content of the underlying sine-Gordon 
field theory in the exponents of the expansion of S a 
on the lattice, whilst these operators do not enter in 
the scaling limit ©. We do not have a satisfactory 
understanding of this result, but wc believe that an 
approach similar to that of RefJ^ might clarify the 
point. 

Since § < no < ir, the exponent of this new cor- 
rection ranges between 1 and 2 and always domi- 
nates over the £ ~ 2 term in ([51]) and competes with 
the correction £~ 2 / Q for a < 2. Notice, however, 
that for mo = — 2/cos^Oi r (wo,A t o) = and thus 
D a (mo, Ho) — 0: this new correction disappears. 
This is precisely the slope that corresponds to the RG 
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flow wc considered before. Also, g'(7r/2) = and thus 
the coefficient in front of this new correction vanishes 
at the crossing point between the free and the bound 
state. Thus, we can conclude that this new term is 
turned on by the existence of a bound state and it is 
a clear signature of a truly interacting theory. More- 
over, the path one chooses to approach criticality se- 
lects a scaling limit in which irrelevant operators can 
be generated and these can modify the perturbative 
series that defines the correlation length, leading to 
something like we observed. 

(c) Straight lines in (7, p) space: Finally, as a general- 
ization of the first case, let us consider the straight 
line 



fj, = fj, + r ■ s , 



(56) 



where r is the slope in the (r, p) plane. Interestingly, 
this trajectory maps into a curve in the ( J z , Jj,)-planc 
approaching the critical point with zero derivative, 
i.e. with a purely quadratic relation in a small enough 
neighborhood of the conformal point (see for instance 
the green line in Fig. @| . Putting (|56|) in (pf5|) we have 



|=4 ff (Mo)^ 2 + 4r£ fl '( M o) li 

We notice the appearance of a strange logarithmic 
correction in the expansion. Inverting (|57|) and plug- 
ging it into the entropy we get 



12a 



ln£ + A a (p) 



E a (r, u) 
ln£ 



(58) 



This kind of logarithmic corrections were found also 
in Rcfiii, but only in studying the limit a — > oo, 
i.e. the so-called single-copy entropy. They were also 
predicted in Ref— as signatures of marginal opera- 
tors in a CFT, but with a different power. From our 
results, it seems that this kind of very unusual correc- 
tions may appear at finite a's, by selecting a proper 
path approaching the critical point. This can be due 
to the presence of a marginally-relevant operator in 
the theory, generated when reaching the critical line 
through a zero-slope curve and corresponding to a 
changing of the compactification radius. We will at- 
tempt a field theoretical analysis of this kind of term 
in a future work. 



V. CONCLUSIONS AND OUTLOOKS 

By using the example of the integrable XYZ chain, we 
proved that, for a massive model, the study of the cor- 
rections to the entanglement entropy as a function of the 
correlation length requires a separate analysis from the 



one that yields the entropy as function of the subsystem 
size. 

For the bipartite Rcnyi entropy of the XYZ model of a 
semi-infinite half-line we found, in the scaling limit and 
as a function of the correlation length, the universal form 
([6|) , where all subleading contributions are explicitly writ- 
ten, thanks to a novel formulation of the reduced density 
matrix in terms of ^-products. We argued that these cor- 
rections are best interpreted in light of a previously un- 
noticed bulk Ising structure of the CTM formulation of 
the model. This means that corrections as a function of 
the correlation length have different exponents compared 
to those depending on the length of the subsystem, un- 
like what was expected from previous studies. This also 
implies that the coefficient of the logarithmic lead- 
ing term has the same value both for c = c = 1/2 of the 
bulk Ising formulation in the mass parameter and for the 
c = 1, c = of the critical chiral free boson model in the 
subsystem size. 

In this respect, it is also interesting to note that the 
reduced density matrix p of (|2Tj) can be written a a 24 i 35 ~ — : 



p <x e 



(59) 



where (r/j, r/j) are (Majorana) fermionic creation and an- 
nihilation operators for single particle states with eigen- 
(57) value 2ej = 2je (note that Hctm is not the Hamilto- 
nian of the subsystem A). This representation strongly 
supports the interpretation that the c = 1 theory is con- 
structed in terms of c = 1/2 (Majorana) characters. 

In Refs i 12 i 16 it was shown that the first correction in 
the Renyi entropy of a critical XXZ chain as a function of 
the subsystem size i goes like £~ 2K / a : where K is the Lut- 
tinger parameter of the model. This fact has also been 
checked in many other critical c = 1 quantum spin chain 
models via DMRG simulations^—. When going to the 
corresponding massive model, assuming that the opera- 
tors responsible for the corrections remain the same, the 
simple scaling prescription^ 1 - would give a term of the 
type ^~ K / a . The results presented here would then indi- 
cate an improbable fixed value for the Luttinger param- 
eter K = 2. This exponent for the massive XXZ chain 
was observed before (take, for instance, Reffii), but its 
nature has not been discussed. Instead, consistently with 
RefJ£, we found in Sec. II VI that the operator responsible 
for this correction is the energy of the underlying bulk 
Ising model. As pointed out in RefJ^, the leading cor- 
rection in I is of the form £~ 2K / a i n the one-interval case 
and £~ K l a for the half-line, whereas for two intervals, in 
the Ising case, the exponent acquires an additional factor 
of 2, which counts the number of twist fields at the edge 
of the interval 3 ^. It would be interesting to perform a cal- 
culation for one interval with two boundary points in our 
case too, to check whether a doubling of the exponent in 
the correlation length would happen in this case as well. 
We also observe that our results stem from the study of 
the formally conformal structure emerging in the Rcnyi 
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entropy of the ground state. It would be of great interest 
to perform a similar analysis for some excited states. 

We also showed that, if one takes into account lattice 
effects, there is a proliferation of new standard and un- 
usual corrections, which in general are path-dependent, 
and can assume the forms: £~ 2/i , £~ 2h ' a , £-(2-/i) ; 
£-2(i+i/a)^ Qr cvcn i/l n (£), where h is the scaling di- 
mension of a relevant operator of the critical bulk the- 
ory. We conjecture that the last logarithmic correction 
is a consequence of some marginally-relevant operator in 
the theory. In addition each of the previous terms ap- 
pears multiplied by a point-dependent and a-dependent 
coefficient that can vanish in even large parametric re- 
gions, preventing the corresponding unusual correction 
from showing up in the scaling limit of the Renyi en- 
tropy. We should remark that the same analysis, carried 
out in terms of the mass-gap instead of the correlation 
length, would generate even different corrections. These 
differences will be analyzed elsewhere for the XYZ chain, 
but are a general feature of lattice models that have to 
be taken into account in numerical analysis. 

In Refii2 it was given a list of possible subleading con- 
tributions to the Renyi entropy for a CFT. From a naive 
scaling argument, one could expect the same kind of 
terms to appear in a massive theory sufficiently close to 
criticality. However, some of the corrections we observed 
do not fit this expectation. This could be due to strictly 
ultraviolet effects that cannot be captured by a QFT or 
to a need to improve the scaling argument. We believe 
that an analysis similar to that carried out for critical sys- 



tems in Refi^ could be successfully applied to a massive 
sine-Gordon theory, by introducing an ultra-violet cut-off 
(the lattice spacing) to regularize divergent integrals and 
extract the corrections from the counterterms^ . This 
type of check would be dual to what we have presented 
here: while in our approach we started from a lattice 
model to infer its universal behavior, in the other, one 
would start from a field theory to understand the origin 
of the correction. In conclusion, we believe that further 
analysis are needed to provide a consistent field theoret- 
ical interpretation of the corrections arising in massive 
models, especially for the relevancy of such problem in 
numerical studies. 
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